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Abstract 

The complicated interactions in presence of disorder lead to a correlated 
randomization of states. The Hamiltonian as a result behaves like a multi- 
parametric random matrix with correlated elements. We show that the eigen- 
value correlations of these matrices can be described by the single parametric 
Brownian ensembles. The analogy helps us to reveal many important features 
of the level-statistics in interacting systems e.g. a critical point behavior dif- 
ferent from that of non-interacting systems, the possibility of extended states 
even in one dimension and a universal formulation of level correlations. 
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I. INTRODUCTION 



The theoretical formulation of physical properties of electronic systems is usually based 
on independent electron approximation. However, experiments on quantum dots as well as 
extended semi-conducting electron heterostructures show that Coulomb interaction between 
electrons is by no means small [1,2]. The influence of Coulomb interaction is particularly 
strong in presence of a disordered environment, with bulk and mesoscopic systems both 
revealing new features [4,7]. Examples include Coulomb Blockade phenomena, low energy 
anomalies in transport and thermodynamics coefficients, non-Fermi liquid effects in effec- 
tively one dimensional structures, various manifestations of the Kondo effects, the fractional 
quantum Hall effects etc. A theory of physical properties including electron electron (e-e) 
interaction is therefore very much required. 

The presence of interactions along with disorder makes a physical system so complex 
that its properties can be analyzed only by a statistical tool. In past, the statistics of 
single particle states of non-interacting systems with external disorder potential has been 
well-modeled by the Wigner-Dyson ensembles, also known as standard Random matrix 
ensembles (RMT) [5,3,7]. The ensembles have also been used successfully to describe the 
statistical properties of the high energy states of many body systems with no disorder 
e.g. nuclei, atoms etc [5,3,7]. The intuition suggests therefore a possibility of random 
matrix modeling of systems when disorder and interactions are both present. Recent studies 
in this direction have led to two types of random matrix models however both impose 
specific conditions on nature and degree of e-e interactions as well as disorder in the system 
[8,9]. In this paper, we consider a random matrix model suitable for generic conditions for 
both disorder and interactions and show that the spectral statistics can be described by 
a mathematical formulation analogous to that of non-interacting disordered systems. The 
parameter governing the localization — > delocalization transition however turns out to be 
different in the two cases. 

The complexity of many body interactions or external disorder potential associates a 
degree of uncertainty in the exact determination of the Hamiltonian. As a result, some 
(or all) elements of the Hamiltonian matrix behave like random variables. The distribution 
properties of the matrix elements are governed by the nature of interactions and disorder; 
their distributions need not be same, may or may not be correlated and some of them can be 
non-random too. For example, a single particle Hamiltonian with a "white noise" disorder 
potential (that is, no e-e interaction or impurity interaction) can be modeled by a random 
matrix with uncorrelated elements. However the presence of local interactions among im- 
purities or electrons or both will result in a correlated randomization of matrix elements of 
the Hamiltonian. Unfortunately not much information for random matrices with correlated 
elements has been available so far. Our objective in this study is to suggest a way to fill 
in this information gap. Here we show that the eigenvalue distributions of various ensem- 
bles, with correlated matrix elements and a multi-parametric probability density, appear 
as non-equilibrium stages of a Brownian type diffusion process. The diffusing eigenvalues 
evolve with respect to a single parameter which is a function of the correlation strengths 
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between various matrix elements. The parameter is therefore related to complexity of the 
system represented by the ensemble and can be termed as "complexity" parameter. The 
solution of the diffusion equation for a given value of complexity parameter gives, therefore, 
the distribution of eigenvalues, and thereby their correlations, for corresponding system. 
As discussed later, the diffusion equation can be solved by using its analogy with the one 
governing the evolution of Brownian ensembles [10]. The analogy also helps in theoretical 
formulation of the level-density correlations of interacting disordered systems. 

The paper is organized as follows. In section II, we derive the diffusion equation for the 
matrix elements of the Hcrmitian operators governing the dynamics in complex systems. 
For a clear exposition of our technique, we first consider the cases with 2"'^ and 3'''^ order 
matrix elements correlations and generalize to n^^ order correlations later on. The diffusion 
equation for the matrix elements is then used to obtain the equation governing the evolution 
of the eigenvalues in section III. To maintain the flow of the discussion, only relevant steps 
are given in the sections II, III; the details can be found in the appendices. As mentioned 
above, the evolution equation for the eigenvalues turns out to be a very well-known equation 
and calculation of the eigenvalue correlations from the equation has been discussed in detail 
many times. We therefore avoid the repetition but give a brief discussion to keep the 
article self-content. This is followed, in section IV, by a discussion of the apphcation of our 
technique to two well-known systems. We conclude in section V by summarizing our main 
results. 

II. SINGLE PARAMETRIC EVOLUTION OF MATRIX ELEMENTS 

Let us consider a complex system represented hy a, N x N Hermitian matrix H with 
Hki — J2s=ii^y''^^ki;s as its matrix-elements. The subscript s to a variable refers to one 
of its components with j3 as their total number. The parameter j3 contains the information 
about underlying symmetry of the system. For example, for systems with time-reversal 
symmetry and integer angular momentum, the Hamiltonian in a generic representation is a 
real-symmetric matrix which gives f3 — 1. The Hamiltonians for system without time reversal 
symmetry are, in general, complex Hermitian which gives (5 — 2. Due to its Hermitian 
nature, the independent real parameters H^^g, which determine all matrix-elements of H , 
are M = N + N{N — l)/9/2 in number. For notational simplification, let us denote them by 

where /i = {kl; s} is a single index which can take a value from 1 — > M. 

The elements of H describe the overlap between various states of the basis in which H 
is represented. A complicated nature of interactions or presence of disorder in the system 
associates a deterministic uncertainty with the matrix elements H^. As the interaction be- 
tween various states is governed by the nature and complexity of the region, the randomness 
associated with the elements can be of various types. In general, the complexity of a 
region can cause multiple interactions between basis states, resulting in correlations between 
matrix elements. In this section, we consider the cases which can be modeled by an ensem- 
ble described by a probability density p{H, h) oc e~^^^^ with function F{H) as a sum over 
various combinations of the matrix elements of H. 
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A. Correlated Gaussian Case 



A Gaussian ensemble of Hermitian matrices with correlated elements can be described 
by a matrix elements distribution 



M M 

p{H,b)^CeM-T.b,,H,,- h,,,,H,,H,,)]^Cp{H,h) 

Ml=l A'l,M2=l 



(1) 



with C as a constant and h as the set of coefficients h^^ and h^^^j_.^. Here the subscripts to 
a coefficient are indicators of the terms present in the product of which it is a coefficient. 
Further, in X^^^^^j' similar pairs are included only once. 

The distribution parameters &2;^i,/i2 the measures of correlations between pairs of the 
matrix elements: < H^j^Hn^ >— ■ In general, different system conditions can give 

rise to different sets of distribution parameters h. A slight perturbation of the system due 
to a change in its parameters perturbs the matrix elements and therefore the probability 
density p{H, b). In the following, we consider a particular flow in the matrix space generated 
by an operator L (describing the diffusion of the matrix elements with a constant drift and 
confined by a quadratic potential). 



E 



d 



d 



2 dH,, 



(2) 



where — 1 + 5^ with 5^ = 1 for /x = {kk;s) and 5^ = for = {kl;s), k ^ I. The 
parameter 7 is arbitrary, giving the freedom to choose the end of the evolution [10]. The 
choice of the above form of L is motivated by our desire to obtain the equation governing 
the eigenvalue-dynamics in a well-known mathematical form. 

The evolution of the probability density p{H, b) , generated by operator L in the matrix 
space, is related to a multi-parametric flow in 6-space. This can be shown as follows. The 
probability density being a function of both H and b, the rates of change of p with variation 
of matrix elements and parameters b can be given as 



dp 
dH, 



Ml 



^Ml ^ '//ii/<2"MlM2-'-'M2 



dp 



db 



^2 

H^iH^^p 



P 



M1M2 

dp 



db 



-H^^.P 



(3) 
(4) 
(5) 



Ml 



where 77 '^'^^^ 



M1M2 



is the frequency of occurrence of the term if^^ in the combination Hf^^H/j^^ (i.e 



Ww2 ~ 2 for /ii = /i2 and 1 for //i p2- With the help of eq.(3) and (4), a drift in the 
matrix space can be written in terms of a drift in parametric space. 
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Ml ^-^/"l 



Ml ^'^Mi Mi,M2 ^'^Mii"2 

where cij^^ = 6^ and c^^^j = ''?m^m2^miM2- Similarly a diffusion in the matrix space can be 
expressed as a combination of drifts in parametric space, 

f>TT2 ~ (*"Ml ~ '^Mli"l)P~'~ ^ '"Ml '^/i^i M2 -^M2 "I" 5^ '^MlM2'^MlM3-^M2-^M3P C^) 

Ml M2 M2,M3 

~ ("^Ml ~ ~ ^ X/ '"A'l'"Ml/"2^7 X/ '^Mli"3'^A'2M3^r (^) 

M2 '^^/^2 M2,M3 '^"m2M3 

A substitution of above equalities in eq.(2) gives us the following, 

Lp^Tp + Cp (9) 

with (7 = J2ij.i [7+ ^(c^i ~ c^iM2)] ^ generator of the dynamics in the parametric 
space, 

^ = X]./'mqI ^ X! ff^T-l^2'K7 ■ (10) 

Ml ^^Mi Mi,M2 ^^MiM2 

Here /^u^ = T^mi ~ Sm2 9ij.2^h2^ij.ih2 /miM2 ~ 7Cmim2 ~ (1/2) S^//3 Qua^iJ-ii^a^i^^na- 

The above equation appears complicated, with many parametric derivatives present on 
its right side. However it is possible to map the multi-parametric flow in the M-dimensional 
6-space to a single parametric drift in another parametric space, say y-space, consisting of 
variables y^, i = 1 — > M where M — [3 + M]M/2. In other words, the generator T of the 
flow in the 6-space can be reduced to a partial derivative with respect to just one y-space 
variable, say yi. 

T(y[6])p^^|„ (11) 

The desired transformation b y required to convert cq.(lO) into the form (11) can 
be obtained as follows. By using the definition ^ = Sfeli ^a^' ^^^^ ^ various b 
parameters, T{b) (eq.(lO)) can be transformed in terms of the derivatives with respect to y, 



where 



V f ^ + V f (13) 

Ml ^'^Ml Ml,M2 ^'^MlW 
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The eq.(12) can be reduced in the desired form of eq.(ll), if the transformation b ^ y 
satisfies following condition: 



A. = 5, 



kl 



for k = 1 



M 



(14) 



The parameters y as a function of b can now be obtained by solving the set of conditions 
(14), 



yi^l2 H / dbf^^^^zl^^^^G + constant 

Ml MliM2 



(15) 



where G = [E^i 4i^//xi + ^St^Ji^w^] ^ and the set of the functions 2:^^' ^ila ^"^^ 

chosen such that the ratio 



A^)Ah -L ^(1) Ah 



MlM2 



Ml 



M1M2 



G 



(16) 



is a complete differential (see appendix B). Note it is possible that the ratio (16) can be 
made an exact differential for different sets z^^^ which will lead to different solutions of y. 
However any two of such solutions for y are different from each other only by a constant 
(appendix B). 

The conditions (14) further imply that parameters y^, k > 1 behave as the constants of 
the dynamics (generated by L) in the matrix space, 



Vk 



= E/ d^Mi4^^ + ^ / d&^i^.4iM2 + constant 



for A; > 1 



(17) 



Ml 



MliM2 ' 



With E^, 4?/mi + E 



^Ml,M2 ^^Ji^J2^MiM2 



for /c > 1. 



The substitution of eq.(ll) in eq.(9), gives the single parametric evolution of the joint 
probability density p{H) in the matrix space 



d 



dH,. 



dp 



dp 
dyi 



(18) 



where yi is given by eq.(15). As the distribution parameters depend on the complexity of 
the system, y^ can be termed as the complexity parameter. 

The parametric space transformation b ^ y maps the probability density p{H, b) to 
p{H,y{b)). As a result, p depends on various parameters yk, k = 1 N. However eq.(18) 
implies that the diffusion, generated by the operator L in the matrix space, is governed by 
yi only; the rest of them, namely yk, k > 1, remain constant during the evolution. Note it is 
always possible to define a transformation from the set 6 ^ y with yk, k > 1 as constants of 
the dynamics generated by L. This can be explained as follows. A matrix element, say Hy, 
describes how a basis state ipi interacts with state ipj through H. This results in dependence 
of the matrix element correlations and, thereby, of the set b, on the basis parameters e.g. 
basis indices. As the basis remains fixed during the evolution, the suitable functions of basis 
parameters can be chosen to play the role oi yk, k > 1. (Note a similar transformation has 
been used to obtain a single parametric evolution of multi-parametric Gaussian ensembles 
of Hermitian matrices; see [10] for details). 
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B. Non-Gaussian Case with third order matrix elements correlations 



Let us consider an ensemble of Hermitian matrices with a 3''^ correlations among its 
matrix elements, described by a probability density 

M M M 

~p{H, h) = Cexp[- KH^, - E b^^,2H,,H,, - ^ b,,^,,,H^,H,,H,,] = Cp{H, b) (19) 

Again the parameters h^j,^^j,^ij,^ are the measures of correlations among three matrix ele- 
ments: < Hnj^Hn^Hn^ >— Q^^°^'~^ ■ Proceeding as in the Gaussian case, we get 



dp 



^2 M2,M3 



(20) 



— — = -Hf,^H^,^H^,^p (21) 



^^/il/i2M3 



The derivatives ^ and remain same as in the Gaussian case (given by eqs.(4,5)), 

V- rr 9p _ ^ _dp_ ^ dp ^ dp 

where c^i^2M3 = v\i'^%^,:fi^il^l2^^3 ^i^h r]^'^^^^^., as the frequency of occurrence of the term 
if^^ in the combination Hn^H^^Hn^. The second order derivative of p with respect to 
can be calculated by differentiating eq.(20) which, on combining with eq.(21), again leads 
to the form Lp = (T + C)p. Here L is still given by eq.(2) however the generator T now 
contains the first order parametric derivatives as well as their products, 

Ml ^"^^1 Ml,^i2 ^"/il/i2 ;il,/i2,iU3 ^'^MH"2/i3 

, V c c ^ (23) 

/il,M2,/i3M4MB P ^I^;i2/i3 ^"M4/iB 

'with /^j = 7C^j — (1/2) X]/i2 fl'At2 [2C/i2CMl/i2 + '^M2At2Atl (^|j2M2/il ~ ■'")]' //^1M2 ~ 7C/X1//2 ~ 
(1/2) X]/i3 fi'/i3 ['^MM3'^A'2M3 2C^3C^i/i2/i3] S-Ild //il/i2M3 ~ 7*^/^1^2/^.3 ~ S/i4 5'/i4 '^/»2A'4 '^A'lM3/'4 • 

Similar to Gaussian case, eq.(23) can also be reduced to a single parametric derivative, 
namely form (11), by using a transformation from the set h to another set y. The operator 
T in this case transforms as 

where 
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The reduction of eq.(24) for Tp to a single derivative ^ will impose following conditions 
on the transformation 6 — > y: 

Ak = NSki (27) 
+ Bki^O (28) 

By solving conditions (27), i/i can be obtained as a function of the parameters b, 
yi = Y.I d^Mi^i^G- + E / d&MiM24i2^ + ^ / d^wM2M3^i2M3^ + constoni (29) 

/il Ml 1/^2 /11,/J.2,M3 

where G = [E^^ ^V/^i + E^i.^a 41U/mim2 + Emim2M3 ^i2M3^iM2M3]"^- Again the choice of 
the functions z'^^^ is so as to make the ratio 



Z^^Ml M ^ ^M1M2^"M/^2 Z^ ^/ilM2M3 MM2M3 

Ml M1M2 ("1^2/1*3 



G (30) 



a complete differential. However the set z^^^ in this case have to satisfy an extra set of 
conditions given by eq.(29). Again as in the Gaussian case, the parameters y^ , j > 1 behave 
as constants of dynamics for this case too. 



I d&Mi^i^ + E / d?^.iM2^ii.2 + E / d?^.iM2^il.2M3 for ^ > 1 (31) 



Ml " Ml>/i2 " Ml,M2M3 

with arbitrarily chosen functions zjj'^ satisfy the constraint X^mi ft^i + I^mi,M2 ^m4M2/("ia'2 + 
Smi,M2 4iLm3^miM2M3 = for A; > 1 as well as the conditions given by eq.(28). 

Using the b ^ y transformation given by eq. (29,31), the dynamics of p{H) with third 
order matrix element correlations can again be described by eq.(18). Note although the 
evolution is governed by same equation, however the complexity parameters are different for 
Gaussian and non-Gaussian cases. 



C. General Case 

A generalized ensemble of Hermitian matrices with correlated elements can be described 
by a matrix elements distribution p{H) — Cp{H) where 

n 

p{H) = n exp 



■ E K{r) 
p(r) 



(32) 
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with C as a normalization constant. Here symbol p{r) refers to a combination of r elements 
chosen from a set of total M — N + N{N — 1)13/2 elements of upper (or lower) diago- 
nal matrix; note the terms present in a given combination need not be all different. The 
11^=1 Hfjtp implies a product over r terms present in the p^^ combination with coefficient 

hp{r) as a measure of their correlation: < ]Y j^iH^^, >— |^^. The X^p(r) is a sum over all 

possible combinations (total (MY) of r elements chosen from a set of total M of them; the 
sum includes only different combinations. (Henceforth subscript p{r) will be written as p 
only unless details required for clarification). 

Using eq.(32), the partial derivatives of p in matrix space and parametric space can be 
given as 



dp _ 
dH,, 



dp 



db 



'p(r) 



r=l p 



P I n ^^^Pi H grr 

j=l J j=l ^^t^ 



MP, 



(33) 
(34) 



As a result, a drift in the matrix space can be related to a drift in the parametric space: 

(36) 



dHu, „ ,s ^ dbr, 

where the term rj^^l^ = Yjm=i ^(/^Pm! counts the frequency of occurrence of the clement 
Hf^ in the p*^ combination of r elements {6{ppj, p) = 1 if pp^ = p and is pp^ 7^ p. Similarly 
a diffusion in the matrix space is related to a nonlinear flow in the parametric space, 



ay 



-lp(r-) 



E 



r'>{n+2~r)p'{r') 



Op' 



dp dp 
dbp+, dbp'+n 



P ^ - E E 

r'<{n+2-r) p'{r') 



dp 



dbpij^p^^2 



Here the notation A-\-B refers to a combination which contains the elements of both A and 
B. Similarly A — B indicates dropping of all the elements of B from A. Further the notation 

p'^ is used to denote a combination (that is, k^^^ power of H^j). 

Again, the matrix space fiow generated by the operator L (cq.(2)) can be related to a 
parametric flow generated by the operator T, Lp = Tp + Cp, where T is now given by 



^p-E^EE 



=1 p{r) 



E E 

_r'>(n+2-r) p'(r') 



h. 



dlogp 



dbp> 



+ fp 



dp 
dbr, 



(37) 



with V = &«^«'E^rf)^l'^] and/p = -Er'y(r') V^JE^^?M^^'^^^^)] + [E^^i^^^?^(^i'^^- 
l)]6^ + 27[E// '7p^'']&p- here q refers to a combination of ri = r' — r + 2 elements, such that 
g(ri) = p'(r') — p{r) + {kVf. Similarly u,u',v refer to combinations u{ro) — p{r) + 1, 
u{ri) = p'{r') + 1 and v{r2) — p{r) + 2 of Tq = r + 1, ri = r + 1 and r2 — r -\-2 elements 
respectively. 
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The desired transformation b ^ y required to convert eq.(37) into form T{y[b])p = ^ 
can be obtained as follows. The substitution of ^^^7 ~ '^^=^^^3'^ eq.(37) with Dkj = 
transforms T in the form of eq. (24) where and B^k are now given as 

p(r) 

Akiy,b)=J2J2fpDkp (38) 

r p{r) 

Bik{y,b)^Yl hppr{x)DipDkp' (39) 

r,p{r) r',p'(r') 

Again, for desired reduction of T to the transformation b ^ y should satisfy conditions 
given by eq. (27,28). The conditions (27) can then be solved to obtain the variables in set y 
as functions of variable in set b (see appendix B) 

yk= J2 J - ^« + + constant (40) 

r,p{r) 

with Gk — Ski Y,r,p{r) fp z^'^^ as arbitrary functions which make J2r,p{r) dbpZ^'^^[l + 

{G — l)Ski\ an exact differential and satisfy the constraint Y,r,r' Sp(r),p'(r') -^p^-^p'^Vpp'^i — 
(the latter is required by the conditions (39)). 

We find, therefore, that, the diffusion of probability density for ensembles of Hermitian 
operators with correlated matrix elements (any order), is governed by a single parameter yi 
with evolution described by eq.(18). The eigenvalue statistics of the above ensembles can 
therefore be studied by an exact diagonalization of eq.(18). 



III. SINGLE PARAMETRIC EVOLUTION OF EIGENVALUES 

The eigenvalue equation for the matrix H can be given as HU = XU with A as the 
diagonal matrix with eigenvalues Aj of as its matrix elements and U as the eigenvector 
matrix (unitary for complex Hermitian case and orthogonal for real-symmetric case). The 
probability density P{E_, y, b) of finding eigenvalues Aj between Ei and + dEi at a given 
Y can then be obtained from the matrix elements distribution, 

N 

P{E, y[b]) - / n S{Ei - Xi)p{H, y[b])dH (41) 
•' 1=1 

Here E_ refers to a diagonal matrix with elements Ei, ..,Eis[. As the Y = t/i-dependence of 
P in eq.(41) enters only through p, a derivative of P with respect to Y can be written as 
follows 

BP r ^ do 

aF = /n/(^.-A.)^di/ (42) 
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The diffusion equation for p, namely eq.(18), can now be used to rewrite eq.(42) as 



dP 
dV 



h + h 



wfiere 



(43) 

(44) 
(45) 



Tlie integral /i can be simplified by using integration by parts 



N Q N 
n=\ 1=1 



pdH 



(46) 
(47) 



A further simplification of the above equation requires a knowledge of the rate of change of 
eigenvalues of H due to a small change in its matrix elements. The rate can be obtained 
by using the eigenvalue equation for matrix H along with unitary (or orthogonal for real- 
symmetric H) nature of its eigenvectors (see appendix A). Using eq.(A2) of the appendix A 
in eq.(47) we get 



d 



dEr, 



{EnP) 



(48) 



The second term can similarly be rewritten as follows 



d 



d 



E 



r. 



^ 2 dH^dH^ 



E 



9t^ d^Xr. 



, 2 dHl 



Now by using eqs.(A3,A5) of the appendix A, I2 can be expressed in terms of eigenvalue 
derivatives of p, 



d'^P d 
' V dE^ ^ V dE„ 



E 



(3P 



E - E 



(51) 



A substitution of /i and I2, given by cqs. (47,51), in eq.(43) leads an equation describing the 
single parametric evolution of the eigenvalues of ensemble p{H), 
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dP 
dV 



E 



n 



d 



'n 



dE, 



d 



'n 



+ E 




'n 



+ 7^ P 



(52) 



The eq.(52) describes the evolution of the eigenvalue density P{E_,y[b]) = P{K,y,y2, ■■Um) 
due to variation of the parameter Y from an arbitrary initial state, say P(£^,y[6o]) = 
P{Ef), Yo, 1/2, ■■Um) occurring at y = Yq- Note here the parameters yj {j > 1) being constants 
of motion, have a same value for both initial ensemble p{Ho,bQ) as well as p{H,b). The 
evolution of the eigenvalues tends to an steady state in limit |^ — or F — > oo. The solution 
of eq.(52) in the limit is a Wigner-Dyson ensemble [7]: P{E) = Ili<j l^i ~ Ejl^e'^'^k^l 
(thus a GOE for /3 = 1 and a GUE for 13 — 2). Note, under certain conditions, the 
steady state solution may also correspond to the eigenvalue distribution of an ensemble 
of anti-symmetric Hermitian matrices. A knowledge of the solution of cq.(52) can now 
help us in determining the n-levcl density correlations -R„, defined as Rn{Ei, E2^ ..En',Y) = 
mzny. I U.f=n+i dEjP{E; Y-Yq). The first order correlation Ri is also known as mean level 
density and its inverse gives the mean level spacing A of the full spectrum (that is, averaged 
of spacings in full length of the spectrum). By a direct integration, eq.(52) can also be used 
to study the evolution of i?„ with changing complexity of the system. 

Equation (52) is applicable for arbitrary values of the coefficients b; it is therefore valid 
for the case of uncorrelated Gaussian ensembles too. The latter have been shown to be good 
models for non- interacting systems [7]. Within random matrix framework, therefore, we 
find that the energy levels of both systems, interacting as well as non-interacting, undergo 
a same diffusion process with changing system parameters. As a result, the level-statistics 
(and related physical properties) in the two cases can be described by the same mathematical 
formulation. However note that, due to different complexity parameters in general, the rate 
of evolution is different in the two cases. 

The advantages of a single parametric formulation of the evolution of eigenvalues is 
manifold and have been discussed in detail in [10] in context of non-interacting systems (or 
uncorrelated ensembles); the similarity of eq.(52) with that of eq.(17) of [10] allows us to use 
the discussion given in sections I.(D),I.(E), II of [10] for the correlated ensembles too (with p 
replaced by E). However, for completeness sake, we briefly review it here again. The eq.(52) 
is similar to the equation governing the evolution of eigenvalues in the Dyson's Brownian 
motion model which was originally introduced by Dyson to describe the eigenvalue-dynamics 
of an ensemble of Hermitian matrices subjected to random perturbation; the non-equilibrium 
states of this model are known as Brownian ensembles (see chapter 8 of [5]). Later on it was 
shown that when an ensemble Hq (fixed or random) is subjected to a random perturbation, 
of strength ^/Y — Yq, by a standard random matrix ensemble V (described by a probabihty 
density e~^^/'^^^'^ ), the resulting ensemble H = Hq + {\^Y — Yq)V behaves like a Brownian 
ensemble (BE) (see chapter 14 of [5], chapter 6 of [6] and [12,10]). Here Ho and V may belong 
to a same symmetry class, with y/Y — Yq governing the parametric eigenvalue dynamics, or 
different symmetry classes with ^^Y — Yq as a parameter for symmetry admixing transitions. 
The statistical properties of BE depend only on the parameter y/Y — Yq besides underlying 
symmetry and many of their correlations are already known [12]. 

As discussed in [12], the mean level density Ri of a BE changes from an initial state to 
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a semi-circular form (typical of Wigner-Dyson ensembles) at the scale of ^{Y — Yq) ~ NAf 
with Ai{E,Y) as the local mean level spacing at energy E; its evolution can therefore 
be described in terms of the parameter (Y — Yq). However the transition of level-density 
correlations to equilibrium, with (Y — Yq) as the evolution parameter, is rapid, discontinuous 
for infinite dimensions of matrices [12]. For small-F and large A^, a smooth crossover can 
be seen in terms of a rescaled parameter A.[E): 

A{E)^j\Y-Yo\/A^ (53) 

The hmits A — > and A — > oo correspond to the level-statistics approaching the initial 
state and Wigner-Dyson ensembles, respectively. As obvious from the definition of A, an 
intermediate state between two limits occurs when the perturbation a/F^^Tq mixes levels 
in a finite energy-range of r local mean level-spacings, r ~ y/Y — Yq/A^^, < r < N. For 
finite size BE, A varies smoothly with changing Y — Yq which results in a continuous family 
of BEs, parameterized by A. However the level-statistics for the large BE (size N — > oo) 
can be divided into three regions: 

(i) initial regime: (Y — YQ)Af'^ — > f{N~^): If the local mean level spacing A; increases 
with size at a rate faster than that of \^Y — Yq, the perturbation will mix fewer number 
of levels as system size increases. The level-statistics therefore approaches its initial state in 
the infinite size limit. 

(ii) WD regime: (Y — Yo)Aj'^ — > f{N): Due to change in with size N being 
slower than that of ^/Y — Yq), even a small change in complexity parameter in this case is 
capable of mixing the levels in an increasing energy range of many local mean Icvcl-spacings. 
This results in an increasing delocalization of eigenfunctions and Wigner-Dyson behavior of 
level-statistics. 

(iii) Critical regime: (Y — YQ)Af'^ — f{N°) — a= a constant: The perturbation 
in this case mixes only a finite (non-zero), fixed number of levels even when the system 
is growing in size. The finite, non-zero A-value in limit N ^ oo therefore gives rise to a 
third statistics, intermediate between initial ensemble and Wigner-Dyson ensemble, which is 
known as the critical Brownian ensemble (CBE). This being the case for arbitrary values of 
a (non-zero and finite), an infinite family of critical BE, characterized by can occur during 
initial ensemble — Wigner-Dyson ensembles. 

The same evolution equations of P for correlated ensembles and BE imply a similarity 
in their eigenvalue distributions for all F-values, under similar initial conditions (that is, 
P{E, Yq) same for both the cases). As a result, one obtains the analogous evolution equations 
for their correlations i?„ too (see eq.(16) of [12]). The mean level density Ri of a correlated 
ensemble can therefore be given by the mean level- density of a BE with same (Y — Yq) 
value and belonging to a same symmetry class (as that of correlated ensemble). Further 
the analogy of evolutions of higher order correlations (n > 1) in the two cases implies (i) 
the discontinuity of transition of Rn for infinite size correlated ensembles, (ii) a smooth 
crossover of i?„ for finite correlated ensembles. The crossover parameter for correlated 
ensembles can again be defined by eq.(54) where now Y — Yq is the complexity parameter 
of the correlated ensemble and A/ as its local mean level spacing. Note, in the case of d- 
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dimensional disordered systems of linear size L, the number of states in a volume of linear 
dimension ( in d-dimensions is n{Q)C'^ with n(0) as the density of states at Fermi level and ( 
as the localization length. Consequently, the typical energy separation between such states 
is ^i{E, Y) = (n(0)C'^)~^. Similarly the mean level spacing of states in the full length of the 
spectrum is A{E, Y) = (n(O)L'^)^^. For disordered systems, the local mean level spacing A/ 
can therefore be expressed in terms of the mean level density Ri as A; = {L/(YR^^ . 

The parameter A, being a function of the distribution parameters of the matrix ele- 
ments, is sensitive to the changes in the system parameters, due to their influence on the 
system-interactions and their uncertainties. Some examples of such system parameters are 
disorder, dimensionality, boundary and topological conditions, system size etc. The presence 
of disorder randomizes the interactions in the system with degree of disorder affecting the 
distribution parameters b and consequently A. The dependence of A on dimensionality and 
boundary conditions can be explained by using a simple example. Consider a N x N lattice 
with a Gaussian site disorder as well as a Gaussian type random interaction between nearest- 
neighbor sites. The lattice Hamiltonian H, in site representation, is a sparse matrix with 
only [Z + 1)N non-zero, independently distributed matrix elements; here Z is the number of 
nearest-neighbors of a site. Consequently only {Z + 1)N 6-parameters (out of N'^) contribute 
to y. As the coordination number Z is different for different dimensions and boundaries of 
the lattice, the Z-dependence of Y results in its dependence on the dimensionality as well 
as boundary conditions of the system. Further the local mean level spacing is also sensitive 
to the dimensionality as well as the boundary conditions. A variation of any of the latter 
parameters can affect both A; and Y and therefore A; (See also [11] where the dependence of 
A on system parameters is explained by considering an example of Anderson Hamiltonian) . 

The size- dependence of A also plays a crucial role in determining the level statistics of 
the correlated ensemble in the crossover regime. In general, both Y — Yq as well as the local 
mean level density are the functions of system-size which results in A^-dependence of A. 
As a consequence, the level statistics in finite systems smoothly approaches one of the two 
end points, namely, A — > or A — > oo, with increasing system size. However, as in BE case, 
the variation of A in infinite correlated ensembles may lead to an abrupt transition, with 
its critical point given by the condition A = size independent. As in the BE case, the size- 
independence of A at the critical point, results in a level-statistics different from the two end 
points. Note if the size-dependence of A^ in a correlated ensemble remains different from 
that of y — Yo under all complexity conditions, the system will never undergo a transition 
in level-statistics. 

IV. EXAMPLES 

In this section, we consider two examples corresponding to 2"^^ and 3^'^ order matrix ele- 
ments correlations and provide the theoretical formulation for 2-point eigenvalue correlations 
for the cases by using the Brownian ensemble analogy. 
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A. Quantum Hall System 



Let H — Hq + V{r) be the single particle Hamiltonian for a disordered quantum Hall 

system with Hq as the kinetic energy of the electrons and V{r) as a space correlated 
disordered potential e.g < V{r)V{r') >= f{r.r') (with < V{r) >= 0). Using the Lan- 
dau states ipnkir) =< r\nk > (the eigenstatcs of Hq) as the basis, H can be written as 
Hnk;n'k' = (^n^n,n'^k,k' + Kfc;n'fe' where Hnk;n'k' =< nk\H\n'k' > and e„ = (n + l/2)huj as the 
eigenvalues of Hq. The interaction between impurities results in a correlation of the matrix 
elements of V and thereby H [13]: 

/ drdrV:,,, {r)i^n,k, i^Wn^k, {^'^n^k, {r')f{r, r') . (54) 

The parameter A can now be determined if mean level spacing A and the real-space cor- 
relations for the potential V are explicitly known. For example, consider the case when mag- 
netic field B becomes much stronger than the disorder potential. The Hamiltonian matrix 
H in this case is divided into various independent blocks (each corresponding to a different 
Landau level) and the statistics of energy states in each Landau Level can be discussed in- 
dependently [13]. For a Gaussian type disorder < V{r)V{r') >= (V'o^/27^cr^)e~l'■""■'l''/^''^ 
the matrix element correlations in the lowest Landau level n — can be given as 
< Voi-fij Vok;oi >=< Vij Vki > where 

with a'^ = {1 + a'^/ll) as a measure of the correlation length of the potential relative to 
the magnetic length = {h/eBY^"^ [13]. Using the notation Ho^-fii = H^i, the distribution 
parameters of the matrix elements of the Hamiltonian in the Landau level n = can be 
given as 

< Hij > = eoSij (56) 

< Hij.g Hkl-s > — ^o^ij^kl+ < ^y;s^W;s > (57) 

with < Vij.sVki;s > as the correlations between different components of the elements of V, 

< Vij,s Vki-s' >= l[< VijVik > +(-1)^-' < VijVki >] 5ss' (58) 

The distribution of the local Hamiltonian for the lowest Landau level can then be represented 
by eq.(l) with parameters b obtained from eqs. (56,57) (here /xi = {ij;s), H2 = {kl;s)); see 
appendix C for an example. A substitution of the 6-parameters in eq.(15) gives us the 
complexity parameter governing the energy level dynamics in the lowest Landau level. As 
shown in the appendix C by a simple case N = 2, the parameters b and, therefore Y, turns 
out to be a function of the disordered potential a, Vq, system length Ly as well as magnetic 
field B and can be varied by changing any one of them. 
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The above discussion is valid for higher order landau Levels too, with complexity param- 
eter still described by eq.(15) however the coefficients h^^^^ are different for different Landau 
levels (see [13] for the matrix element correlations of potential V). The rate of transition 
of level-statistics therefore differs, in general, from one Landau level to another. For weak 
magnetic fields, where various Landau levels can not be considered as independent, H can 
still be represented by the ensemble cq.(l). However, now the number of coefficients h^^^^ 
which contribute to Y is much larger (due to correlations between levels in two different 
Landau Levels). 

In absence of disorder, under the independent Landau Level approximation, all energy 
levels in a given Landau Level are degenerate and matrix H = Hq is a diagonal matrix with 
a Poisson behavior for its eigenvalues (due to dominance of zero spacings). The switching 
of disorder removes the degeneracy and delocalizes the wavefunction if the impurities are 
interacting. The degree of delocalization depends on the strength of impurity interactions 
with respect to the magnetic field strength B. If the latter is strong enough to mix the 
levels in an energy range of many mean level spacings, (which corresponds to the limit 
H !^ V), the energy levels of H show a GUE behavior. (This is similar to the case of 
strongly interacting many body systems e.g. the statistics of resonances in complex nuclei 
which can be well- modeled by GOE or GUE [7]). Under an intermediate state of disorder, 
therefore, the ensemble H lies between the Poisson ensemble and GUE and can be modeled 
by eq.(l). The level statistics for this case can then be given by the one for a BE appearing 
during a Poisson — > GUE transition . The two point correlator R2 [12] for states in lowest 
Landau level can therefore be given as 



where R2{r, 00) = 1 — '^'"2 2'' (the GUE limit); the above formulation was obtained for the 
BE in [12]. 



Let us now consider the case with weak magnetic fields where the interaction between 
various Landau levels can not be ignored. The eigenvalue spectrum of Hq for this case 
behaves as a uniform spectrum [13] (that is, an initial spectrum of uniform spacing). The 
switching of interacting-impurities potential V again results in broadening of the wavefunc- 
tions. In the limit where impurity interactions are strong enough to mix energy levels in 
different Landau levels (that is H V), the eigenvalues of H show a GUE behavior. The 
varying degree of the interaction between impurities therefore leads to a transition of the 
level statistics from uniform spectrum GUE behavior. The two point correlation for the 
level-statistics at any intermediate stage of impurity interaction can then be given by that 
of a BE appearing during the uniform spectrum — > GUE transition: 




dx / dz cos(27rra;) exp — 87r^Aa;(l -|- a; -|- 1z\fx) 





q=—oo 



with A ^ 00 corresponding to GUE limit. 
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B. Disordered Systems with fermionic interactions 



Consider a general Hamiltonian for spinless interacting fermions 



(61) 



ij ijkl 



here the states \i >= af\0 > describe a fixed basis of m single-particle states with Vij 
as matrix elements of the one body Hamiltonian and Uijki as the antisymmctrized matrix 
elements of the two body interaction U. The presence of disorder randomizes both V and 
U. For example, for F as a white noise, its matrix elements are independently distributed 
random variables. However the fermionic interaction results in correlations among matrix 
elements of U : 



< UijkiUkimnUmnij >= / d'^rid'^r2d'^r3'ip*j{ri)i!ki{ri)ijli{r2)ijrnn^ I 



with /(ri, as the interaction between two fermions at positions ri and r2. The Hamilto- 
nian H will therefore be a matrix with varying degree of correlations between its elements 
and can be represented by ensemble eq.(19). In absence of fermionic interaction, the matrix 
H = V and its statistical properties depend on the degree of disorder. For example, in 
presence of strong disorder, the eigenvalues of V show a Poisson distribution with localized 
eigenfunctions [11,3,7]. The weak disorder hmit of V shows a Wigner-Dyson distribution for 
its eigenvalues with extended eigenfunctions. Similarly in absence of disorder, H = U and 
its statistical behavior is governed by the fermionic density in various parts of the system. 
The presence of almost uniform fermionic density in the system leads to delocalization of 
eigenfunctions and U behaves like an ensemble of anti-symmetric Hermitian matrices (see 
[5] for details on anti-symmetric random matrices) . However non- uniform electronic density 
in various parts of the system (that is, stronger interactions in certain parts of the system 
as compared to others) can result in localization of wavefunctions and thereby a Poisson 
behavior for the eigenvalues of U. 

In presence of both disorder as well as fermionic interaction, the behavior of H is gov- 
erned by the inter-competition between them. In this case, it is preferable to represent 
H in the = rri^ dimensional basis of two particle states \ij >; the choice of the basis 
results in appearance of the disorder elements as the diagonal elements of H and fermionic 
interaction elements as the off- diagonals: < ij\H\ij >= Vij and < ij\H\kl >= Uijki- The 
changing strength of the fermionic interaction in presence of disordered potential subjects 
the level-statistics to undergo a transition from the initial state (given by the statistics of 
V) to Wigner-Dyson statistics (when U c:^ V). However, in the limit when the disorder 
potential becomes negligible as compared to fermionic interactions, the level-statistics of H 
approaches that of anti-symmetric matrices. As obvious, the level-statistics for all other 
cases, corresponding to different strengths of disorder potential and fermionic density, will 
lie on the transition curve from Poisson — > Wigner-Dyson — > anti-symmetric ensembles. 
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Using eq.(19) as the model for intermediate states, the level-statistics for this case can then 
be described by BE lying between Poison — > anti-symmetric ensemble. 

The ground state properties of electrons in nano-particles or quantum dots, that is, finite 
systems of fermions interacting via coulomb forces, are not yet fully understood. Using 
ensemble (19) or (32) as their model, the physical properties of these system can now be 
probed further and earlier experimental observations may be explained. For example, it has 
been experimental observed that the peak-spacing statistics for an irregular quantum dot in 
Coulomb blockade regime undergoes a crossover from Wigner-Dyson to Gaussian behavior 
as the strength of electron-electron interaction increases [1,2]. Within our formulation, 
the observed behavior can be explained as follows. The single electron dynamics inside a 
quantum dot of irregular shape is chaotic; the level-statistics of single particle Hamiltonian 
V can therefore be modeled by the Wigner-Dyson behaviour [3]. The addition of more 
electrons switches the potential U, however, due to non-uniform electronic density during 
initial stages, the correlation between various matrix elements of U need not be same. The 
statistical behavior of the quantum dot can therefore be described by an ensemble given by 
eq.(19). As electron density increases, U dominates over V and the level-statistics approaches 
the behavior of anti-symmetric matrices. It is already known that the nearest-neighbor 
spacing distribution for non-central spacings in the spectrum of anti-symmetric Hermitian 
matrices behaves like a Gaussian distribution [5]. The observed Gaussian behavior of the 
peak-spacings in the strong interaction limit is therefore well in agreement with theoretical 
expectations. Using ensemble (19) as a model for the quantum dot Hamiltonian, the behavior 
of peak-spacings in the intermediate regime can be predicted to be similar to that of a 
BE appearing during a cross-over from GOE — > anti-symmetric ensembles. A detailed 
quantitative analysis of such cases is still in progress and will be published elsewhere. 

In general, the size dependence of the parameter Y — Yq and the local mean level spac- 
ing for a system with e-e interaction is different from the non-interacting systems. 
In interacting systems, therefore, the critical point of level-statistics, given by condition 
A = — independent, can occur at a disorder strength (or energy) different from the one 
for non-interacting systems. This is consistent with the results given by renormalization 
group techniques [15] which show that the introduction of interactions into quantum dots 
can produce phase transitions in the limit of weak disorder, leading to behavior qualitatively 
different from the non-interacting case. Further for one dimensional non-interacting systems, 
it is known that the wavefunctions are localized even in a weak disorder limit. However our 
formulation indicates that the fermionic interaction may lead to extended states even in one 
dimensional disordered system; the implication is in agreement with earlier studies in this 
context [16]. 

V. CONCLUSION 

In summary we show, for the first time, that the level-statistics of disordered systems with 
interactions is governed by a single parameter, namely, the rcscalcd complexity parameter A. 
Note the level-statistics of non-interacting systems can also be described in terms of A [10]. 
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However the introduction of interactions modifies the dependence of A on system parameters 
which can significantly affect the location of the critical point for the phase transitions and 
corresponding level statistics. Our study also reveals a deep level of universality underlying 
physical systems, namely, the Brownian ensembles as the statistical back bone of both 
interacting as well non-interacting systems. This universality should be explored in full 
detail as it may reveal many new connections among a wide range of complex systems and 
can be helpful in theoretical formulation of many of their physical properties. 



APPENDIX A: THE CHANGE OF EIGENVALUES AND EIGENFUNCTIONS 



The eigenvalue equation of a complex Hcrmitian matrix H is given by HU — UA with 
A as the matrix of eigenvalues A„ and U as the eigenvector matrix, unitary in nature. As 
obvious, a slight variation of the matrix elements of H will, in general, lead to variation 
of both the eigenvalues as well as the eigenvectors and associated rates of change can be 
obtained as follows (see appendices A-E of [10] for more details): 

As A„ = J2i.j UniHijU*j, the rate of change of A„ with respect to H^i-s (with s referring 
to real, s = 1, and imaginary, s = 2, parts of Hki) can be given 

^ = ^[UlnUL - {-imnUU (Al) 

Orlklis Qkl 

where gki — '^-\- Ski- This can further be used to obtain the following relations 

2 a\ 

E E nn^H^i,s = E HkiUinUl^ = Xn (A2) 

k<ls=l'^^kl;s k,l 

and 

2 

E 9kl E ^77^ -^rf- = 2(^mn (A3) 
k<l s=l "^kl;s OHkl-s 

As obvious from eq.(Al), the second order change of an eigenvalue with respect to a ma- 
trix element requires a knowledge of the rate of change of one of the eigenvector components 
with respect to Hki. The latter can again be obtained by using the eigenvalue equation, 

^ = ^ E Y^^^-^^*^rnUln + {-ir'U;^Ukn) (A4) 
Onkl;s 9kl jn^n ~ 

Now by differentiating eq.(Al) with respect to a^^s a-nd by using the eq.(A4) we can 
show that 

E^.^ES^ = 2/3E^r4^ (A5) 

k<l s=\ '-'^^kl;s m 

For the real-symmetric case, the corresponding relations can be obtained by using C/+ = 
U'^ (as eigenvector matrix is now orthogonal) in eq.(Al) and taking Hij.2 — for all values 
of 
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APPENDIX B: SOLUTION OF EQ.(38) 

According to theory of partial differential equations (PDE) [14], the general solution of 
a linear PDE 

^Pi{xi,X2, ..,xm)-k— ^ R (Bl) 

1=1 * 

is F{ui,U2, ..jUn) = where F is an arbitrary function and Ui{xi,X2, ..,Xn, Z) = Ci (a 
constant), i = 1, 2, .., n are independent solutions of the following equation 



dxi dx2 dxfc dxM dZ 



(B2) 



Pi P2 Pk Pm R 

Note the function F being arbitrary, it can also be chosen as 

F ^ Y^^u, - c,) = (B3) 

The equations for various yj in the set of eq.(38) are of the same form as eq.(B2) and, 
therefore, can be solved as described above. Let us first consider the equation for 7/1; its 
general solution can be given by a relation F{ui, U2, .., um) — where function F is arbitrary 
and Uj are the functions of M parameters of set b such that Uj{{b},yi) = Cj (with c^'s as 
constants) . The functions Uj are the independent solutions of the equation 

fp(2) fp(r) 

where the equality between ratios is implied for all possible combinations p{r) of r terms, 
r = 1 — > n, with M as the total number of combinations. It is easy to see that each of 



the above ratios is equal to ^^^^ " -where z^^^ are arbitrary functions. The eq. (B4) can 

therefore be rewritten as 



■ir,p 



dyi = T ' g), ^ (B5) 

I Zn tin 



A solution, say Ui of eq.(B5), or alternatively eq.(B4), can now be obtained by choosing the 
functions z^^ such that the right side of the above equation becomes an exact differential: 

Ui = yi— ^ j dbpZ^^Gi — constant (B6) 

r,p(r) 

where Gi = [I]r,p(r) ^p^"*^?]^^- The general sohition for yi can therefore be given by a 
combination of all possible functions u obtained by using arbitrary set of 2;-functions. It 
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can be shown that each such solution differ from the other only by a constant: uj = Ui + 
constant; (this is due to equahty of the two ratios obtained by choosing two different sets 
2;^^) of the functions). The yi can therefore be written as follows, 

yi^ J dbpZ^^^Gi + constant (B7) 

r,p{r) 

which gives eq.(40) for k = 1. 

The set of equations (38) can similarly be solved for other yj {j > 2). For example, the 
solution of eq.(38) for y^ can be given by the function Fk{vi, ..vm) = where Vj{{b},yk) = 
constant are the independent solutions of following equality 



d&p(i) _ _ d6p(2) _ dbp^r) _ dy, 



(B8) 



fp(l) fp{2) fp(r) 

A solution, say vi, of eq.(B8) can now be given as 

vi = yk — 'Y j dbpZ^^ — constant (B9) 

r,p{r) 

where z^'''> are arbitrarily chosen M functions which satisfy the condition 

E = (BIO) 

r,p{r) 

As obvious, one possible choice for z^''^ functions satisfying the above condition is z^'''^ — 
for all p{r) which gives yk = constant. 

As each solution of eq.(B8) is different from the other only by a constant, the general 
solution for yfe, k > 1, can now be given as 



^ J dbpZ^J"^ + constant (Bll) 

r,p{r) 

The eq.(B7) and eq.(Bll) together give the set of eqs.(40). 



Vk = 

r,p{r) ' 



APPENDIX C: EXAMPLE FOR QUANTUM HALL CASE 

Within independent Landau level approximation, the Quantum Hall ensemble can be 
described by the probability density p{H) given by eq.(l) with coefficients b given by 
eqs. (57,58). However it can further be simplified by choosing the origin of energy at eg 
which makes b^ =< >= 0. The matrix clement distribution in the Quantum Hall case 
can now be described by a probability density 

M 

p(i7, 6) = Cexp[- b,,,.H,,H,.,)] (CI) 
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with C as the normahzation constant 



= / 

Here B^^^ is the matrix of coefficients h 



correlation < 11^^11^^ >: 



n Det5,)-V2 

s=l 

^,/,2- The parameters ^2 



(C2) 

are related to second order 



1M2 



91nC 



2 96 



(C3) 

(C4) 



where r7y:iijU2 = 2 and 1 for pairs {//i} = {//2} and {//i} ^ {112} respectively. Let Q^^") be 
the matrix with its elements as the correlations between elements of if, that is, {Q^^^) ij,iti2 — 
(4/^wM2) < H^^H^^ >. By using eq.(C2) in eq.(C4), we get 



91nDet[5(^)] Cof(5 



96^,^2 Det[S(-)] 



-1 



IJ'2lJ-l 



(C5) 



where Coi{B'^^^^) implies the cofactor of the element B^-^^^ in the matrix B^. This implies 
Q{s) ^ The matrix B^'^ for the QH case can therefore be obtained by inverting the 

correlation matrix Q^"*-* = {Q^/^in^}- 

Let us consider the case N — 2. Using the notation, various components of matrix 
elements can now be denoted as Hi = Hu-i, H2 = ifi2;i, H3 = if22;i, H4 — ifi2,2- Following 
eq.(60), the correlation matrix Q^^^ in this case is a 3 x 3 matrix 



/ a 



o 2ax \ 



o a/2xi 



(C6) 



\2ax 



a I 



where a = {y^l2lJjyOL^/2^), x = e'^/^"" and xx = e^^'/^. By using the relation (C5), S(^) 
can be given as 



a(l - 4x2) 



/I o -2x \ 

o 8x1(1 - 4x2) 



\-2x 







(C7) 



Due to Hermitian nature of if, only its off-diagonal elements have imaginary parts. For 
N — 2 case, therefore, S^^) is just a 1 x 1 matrix, corresponding to correlation < H/^H/i > 
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with = (12; 2): B^'^^ = 644 = (2 < H4H4 >) ^. The parameter set y for this case can 
now be obtained by solving the condition A/. = Ylij=i fij^^ = ^fci where fij = "ybij — 

(1/2) J2k=i 9kCikCjk', here Qk — '^ for odd k and Qk — ^ for even k. As discussed in appendix 
(B), a solution of the condition can be given as 



d&ii _ d6i2 _ d6i3 _ d&22 _ d^si _ d623 _ d&33 _ d644 _ dj^ ^^^^ 

/ll /12 /l3 /22 /3I fz2 /sS /44 ^fcl 

The above equations can now be solved by making the ratio exact differential, (following 
from eq.(B5) of appendix B) 

^ = ^ (C9) 

Gk Okl 

where dF^ = Z^jj 4j ""dfeij and Gk = J2i,j ^ij '' fij with z^'^^ as arbitrary functions. Using 
eq.(C7), the Gk can be shown to be 

Gk = _ [(4? + 4?)[(7a - 4) - 4x2(7a + 1)] + {z^ + 4?)4x(8 - 7a(l - 4x2)) + +16zg^a-'x^ 

Similarly 

dFk = E 4?<ibij = [(4i^ + 4?)dx2 - 2(4? + 4?)d(xx2) + 2(44? - 4?)d(xi/a)] (Cll) 

where X2 — a~^{l — 4:x'^)~^. As both dFi and Gi are functions of Ic and a (through a and 
a), yi will turn out to be a function of parameter Ic and a, 



- Gi J Gi{lc,a) 
Proceeding similarly for /c > 2, a solution for y^ can be given as 



yk^j:f d6.,4f (C13) 

where z'^^^ satisfy the conditions = 

Note the condition G2 = is satisfied for a following choice of z'^^^ : z^^x = ~43^ ) -^la"* = 
—4? and 42'* = 41'* = 0- The 2/2 for this choice turns out to be a constant. Similarly the 
condition G3 = can be satisfied for a following choice of z^^^: ^]k^^fz\ = ^'^^.^^"■'■/"t^,-,^ 2 

and ^ = " ^l{^ax^^-i%)^ ■ Using these z values in eq.(cl3), one can obtain y^ as a function 
of a and Ic- Note although y^ varies with changing a and Ic-, however Y^ij fij§^ = and 
therefore y^ does not affect the evolution of p{H). 
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